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The concept of "age" as a parameter for the description of the state of development of high 
energy showers in the atmosphere has been in use in cosmic ray studies for several decades. In this 
work we briefly discuss how this concept, originally introduced to describe the average behavior 
of electromagnetic cascades, can be fruitfully applied to describe individual showers generated by 
primary particles of different nature, including protons, nuclei and neutrinos. Showers with the same 
age share three different important properties: (i) their electron size has the same fractional rate of 
change with increasing depth, (ii) the bulk of the electrons and photons in the shower (excluding high 
energy particles) have energy spectra with shapes and relative normalization uniquely determined 
by the age parameter, (iii) the electrons and photons in the shower have also the same angular and 
lateral distributions sufficiently far from the shower axis. In this work we discuss how the properties 
associated with the shower age can be understood with simple arguments, and how the shapes of 
the electron and photon spectra and the relative normalization that correspond to a certain age can 
be calculated analytically. 

<D • PACS numbers: 96.50.S-, 96.50.sd, 13.85.Tp 
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^ i I. INTRODUCTION 

Q . The concept of the "age" of a shower has been in use in the cosmic ray community for more than half a century. 
^ ' The concept first emerged [ij in the study of the average longitudinal development of purely electromagnetic showers 
c/2 , generated by photons or electrons. It was then also applied |3j to the lateral distribution of electrons around the 
shower axis. Soon, it was also understood that it is possible and useful to assign an "age" also to individual showers, 
and that the concept is applicable also to showers generated by hadronic primary particles such as protons or nuclei. 
—i • Some recents works have rediscussed the concept of of shower age for the showers generated by ultra high 
K*" I energy cosmic rays in the Earth's atmosphere. Giller et al [H and Nerling et al. [H] have studied with montecarlo 
. methods the showers generated by high energy protons and nuclei in air, and have observed that the energy and angle 
distributions of the electrons (in this work with "electrons" we will refer to the sum of electrons and positrons) in the 
showers have shapes that to a good approximation are only determined by an age parameter s defined as: 



O' S(i, tmax) = T—r- ■ (1) 

OO ' 
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where t is the depth in unit of radiation lengths and tmax is the depth where the shower reaches its maximum size. 
The shower size is defined as the total number of charged particle integrated over all energy, and eS'ectively coincides 
ij - with the electron size. These results have been extended to the lateral distribution of electrons by Gora et al. 
rS ' This property of "universality" is obviously very important for the analysis and interpretation of high energy cosmic 

ray observations, and it is therefore desirable to have a deeper understanding of its origin and of its limitations. 
■ ■ In this work we want to review critically the concepts of shower "age" and "universality" . One of the main points 

we want to make is to argue that the definition of age of equation ([T]), while reasonably accurate in most cases, is in 
general not correct, and should be replaced by a better motivated and more accurate definition. 

The essence of the concept of shower age can be understood observing that all showers at the maximum of their 
development are in an appropriate sense "similar" to each other (that is have the same "age"). This "similarity" is 
represented by the fact that in all showers at maximum the energy spectra of "most" electrons and photons have 
the same shape and the same relative normalization. These particles also have the same angular distributions (that 
are obviously strongly correlated with energy) and, for an equal density profile of the medium where the shower is 
propagating, also the same lateral distribution around the shower axis. 

The idea that all showers at maximum, that is at the depth where the derivative of the shower size N{t) vanishes, 
are "similar" independently from the energy and nature of the primary particle, can be naturally be generalized, 
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stating that all showers that have the same fractional rate of change with depth, that is the same "size slope" A: 

A- ^ ^^^'^ (2) 

are also "similar" . This means that to each value of the A correspond well determined shapes of the electron and photon 
spectra (again only valid for "most" particles) and a well determined relative normalization for the two populations. 
It is intuitive (and will be later verified by detailed calculations) that the energy spectra of electrons and photons 
become progressively softer as A decreases going from positive (when the shower size grows) to negative (when the 
shower size decreases) values. 

There is a one to one mapping between the values of the size slope A and the values of the shower "age" s. This 
mapping is encoded by a function that, in the notation introduced by Rossi and Greisen in their "classic" paper 
is called Ai(s). The general definition of the shower age is therefore: 



.Nit) dt J 



(3) 



where Aj~^ is the inverse function of Ai(s). The function Ai(s) (that will be discussed in more detail in the following) 
is monotonically decreasing and has a single zero at s = 1 , therefore according to equation ([3]) showers have age s — 1 
at maximum and age s < 1 (s > 1) before (after) maximum. 

The definition of the age parameter ([3]) may seem at first sight (and in some sense actually is) arbitrary, since the 
size slope A itself or any monotonic function of A are also perfectly adequate to identify "similar" showers. The choice 
of the particular mapping of equation ([3]) is motivated by the fact that one can attach a direct physical meaning to 
the quantity s. The shapes of the electron and photon spectra for E above the critical energy e (the electron critical 
energy e is the average energy lost by an electron in a radiation length, and corresponds also to the electron energy for 
which radiative and collision losses are equal) and E ^ Eq (with Eq the primary particle energy) are well represented 
by a power law: 

n,iE) ^ n^{E) ^ E-^'+'^ (4) 

These power law behaviors stops when E approaches (from above) the electron critical energy. For energies below 
iJ ~ e the electron spectrum has a sharp cutoff, while the photon spectrum has a "knee" and takes the form E^^. 
The precise shapes of the cutoff of the electron spectrum, of the knee of the photon spectrum (that is the transition 
from the form to the E^^) and the relative normalizations of the photon and electron spectra are all entirely 

determined by the age s (or equivalently by the size slope A) and can be computed in detail. 

The commonly used definitions of age in equation ([T]) always coincides with the general definition ^ at shower 
maximum and therefore, by construction, it is a reasonably good approximation for showers sufficiently close to 
maximum. The motivation for the more general definition may appear as only a formal question of "principle" . In 
fact in some circumstances the two definitions are significantly different, and the general definition gives the correct 
shower age. The definition ([1]) coincides with the correct one only for a particular shape of the shower longitudinal 
that is known as the "Greisen profile" Q . In fact Greisen profile and the age definition ([T]) are intimately connected, 
and can be seen (with the mediation of the function Ai) as the integral and the derivative of each other. The Greisen 
profile (discussed below in section HV]) describes accurately the average development of purely electromagnetic showers, 
but is only a rough approximation for the description of individual hadronic showers, and naturally fails completely 
in the description of neutrino-induced showers. The deviations of the definition ^ from the true age ^ are of the 
same order of the deviations of the profile of a shower from the Greisen profile that has the same imax- 

The authors in U Q have calculated with montecarlo methods the shape of the electron spectra in hadronic showers 
of different age. The parametrizations of their results are essentially identical to the shapes of the electron spectra in 
showers of the same age calculated several decades ago by Rossi and Greisen [l| . These modern works have therefore 
effectively only "rediscovered" with montecarlo methods what should be called the "Rossi-Greisen" spectra. We want 
to attract attention to this fact for three reasons. The first one is that it is obviously appropriate to give credit to 
the remarkable work of the pioneers. The second is that the works of [1, Q do not include a discussion of photon 
spectra. The shapes of these spectra and their relative normalization with respect to the electron ones are also 
determined unambiguously by the shower age, and have been also computed explicitely by Rossi and Greisen. Finally 
the derivation of the spectral shapes obtained by Rossi and Greisen with analytic methods allows physical insights 
on the origin and limitations of the "universality" of the spectra, that are not easily deducible from a montecarlo 
calculation. 

It should be stressed that the "universality" of properties for cascades of the same age has clearly limitations since 
it only applies to "most" but not all particles in the shower. For example, the "similarity" among showers at the 
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maximum of their development, (that is at age s = 1) does not imply that the showers simply differ by the absolute 
normalization of their electromagnetic component. Considering at first the case of purely electromagnetic cascades, 
at maximum the showers generated by a photon of initial energy Eq contain (in essentially all cases) more high energy 
particles than the showers generated by photons of lower energy. These high energy particles are negligible in number 
and do not contribute significantly to the total size but in general carry an important fraction of the shower energy, 
they "feed" the shower development and influence its development. The showers generated by other types of primary 
particles have "cores" of different structure and particle content, and follow different development profiles. 

This work is organized as follows: in the next two sections we review a very well known subject, discussing the 
average longitudinal evolution of purely electromagnetic showers first in "approximation A" , that is neglecting the 
electron ionization losses, and then in "approximation B" . The concept of age emerged naturally in these studies. 
In approximation B the shower equations have "elementary solutions" labeled by the parameter s, these solutions 
correspond to the "universal spectra" of showers with age s. The following section discusses the well known "Greisen 
profile" that describes the average longitudinal development of purely electromagnetic showers. Finally we discuss 
the evolution of individual hadronic showers, and give some conclusions. 



II. ELECTROMAGNETIC SHOWERS IN APPROXIMATION A 



The evolution of purely electromagnetic showers can be studied [l| using two sets of simplifying assumptions called 
"Approximation A" and "Approximation B" . 

In approximation A the only processes considered for the shower development are pair production for photons, 
and bremsstrahlung for electrons. The differential cross sections for these processes are described by the asymptotic 
formulae valid at high energy. The electron energy losses due to collisions with electrons and nuclei of the medium 
are neglected. 

The average longitudinal development of electromagnetic showers is described by the two functions ne{E,t) and 
n.y{E,t) that give the differential energy spectra of electrons and photons at depth t. In this work we are following 
the notation introduced by Rossi and Greisen in however here we introduce different symbols. Rossi and Greisen 
indicate the differential (integral) electron spectrum as Tr{E,t) {n{E,t)) and the photon spectrum as j{E,t); the 
subscript notation used here is more suitable to the extension of the formalism to hadronic showers where other 
particle types are present. 

In approximation A the evolution of the electron and photon differential spectra is described by the two integro- 
differential equations: 



dne{E,t) 
dt 



dv ipo{v) 



ne{E,t) - 



1 - V 



lie 



E 
1-v 



^ du 



E 



(5) 



dnJE,t) ['^ dv , , f E \ ,^ , 

' = — ^{v)nJ-,t\ -aon^{E,t) . (6) 

In the right hand side of equation ^ the first term describes the (e — > e) contribution, and the second one the (7 e) 
processes. In the right hand side of equation ^ the first term describes the (e — s- 7) contribution, and the second one 
photon absorption. The differential cross sections for bremsstrahlung ip{v) and pair production ipiu), and the photon 
absorption cross section ctq are given in appendix A. 



A. Elementary solutions 

In the system of equations ([5]) and ([6]) no energy scale is present. Accordingly these equations have a set of 
"elementary" , scale invariant solutions of form: 

that are power laws in energy, and change exponentially with the depth t. Inserting these solutions in the shower 
equations ([5]) and ^ one obtains a quadratic equation for A(s) that has the two solutions: 



Ai,2(s) = -I (Ais) + ao) ± \^{A{s) - aof +AB{s) C(s). (8) 
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To each solution corresponds a photon/electron ratio: 

^^''^(^) = — rrVr (9) 



The auxihary functions A{s), B{s) and C(s) appearing in the definitions ([8]) and ([9]) are given exphcitely in appendix 

A. The functions Ai,2(s) are shown in fig. [11 the function r^^'' is shown in fig. El 

Even if Ai(s) is already given by an explicit analytic expression in ([8]), it is very useful to follow Greisen and 
introduce the simpler expression 

Ai(s) = i(s-l-3 1ns) (10) 

that is a very good approximation for Ai(s) with deviations smaller than 2% in the interval 0.6 < s < 1.4. A 
comparison of the exact and approximate expressions for Ai(s) is shown in the top panel of fig. [H The usefulness of 
this simpler functional form will be clear in the following. 

The existence of two solutions Ai.2(s) for each s value is physically simple to understand. If one starts at i = 
with populations of electrons that have power law form with the same slope, but arbitrary normalizations, the spectra 
mantain identical power law shapes at all t, but change their relative and absolute normalization. The spectra reach 
first an asymptotic 7/e ratio with a t scale |A2(s)|~^, and then evolve exponentially oc e'^^^^^* mantaining a constant 
equilibrium ratio. The convergence to an asymptotic 7/e ratio is the fastest process since |A2(s)| > |Ai(s)| for all s 
values. 

As an explicit example, an initial power law, pure electron spectrum: 

f neiE,0) = K 



evolves in t as: 



{ n^{E,0) = 

n,iE,t) = [(Ai(.)+ao)e^^(^)*-(A2(.)+ao)e^^(^)*] E-(^+^^ 



(11) 



(12) 



The spectra remain power laws for all t values. For t ^ |A2(s)|^^ one can set e'''^^^)* to zero, and the spectra evolve 
in t as a simple exponential (oc e'*'^^'*^*) with an asymptotic 7/e ratio C(s)/(Ai(s) + ctq) that corresponds to the first 
solution in ([9]). 

In summary, the solutions 

neiE,t) = K' e^i(^)* 

(13) 

n^{E,0) = K' e^i(^)* r\^\s) 

is a sort of "attractor" , and any combinations of photon and electron spectra power laws spectra of the same slope s 
"converge" to this solution. 

The only t independent solution corresponds to s = 1 and is particularly important: The existence of this t- 
independent solution: 

ne{E,t) = KE-^ 

n^{E,t) = Kr^ E-^ ^ ' 

and its energy dependence oc E~^ do not depend on the detailed form of the pair production and bremsstrahlung cross 
sections, and can be immediately understood observing that a power law spectrum of form E~'^ contains equal amount 
of energy in each energy decade and that the bremsstrahlung and pair production processes that "mix" the electron 
and photon populations, conserve energy and are scale invariant. Only the value of the 7/e ratio that corresponds to 
this solution depends on the detailed form of the cross sections and is: 

= r^Hl) . . ^ . (1±^ . 1.31 (15) 

(To o-o (7/9 - b/6) 

The fact that Ai(s) is positive (growing solution) for s < 1 and negative (decreasing solution) for s > 1 is also 
independent from the detailed form of the cross sections, and is a simple consequence of the fact that in a power 
law spectrum of form £'~(^+^) the energy contained in each decade increases with E when s < 1 and decreases when 
s > 1. 

This very elementary discussion already illustrates how the i-dependence of the shower development is intimately 
related to the shape of the energy spectra of the particles in the shower. 
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B. Showers generated by a primary 7 or e^. 



The evolution of the shower generated by a primary electron or photon of initial energy i?o, in approximation A 
has been calculated by Rossi and Greisen [l| . Before describing these solutions explicitely we can observe that some 
important properties can be readily deduced with simple considerations. The differential spectra of electrons and 
photons in the solutions have the scaling form: 



na{Eo,E,t) = fa 
Eo 



E 



(16) 



(the subscript a runs over the 4 cases: e — > e 
the scaling form: 

NaiEo, E^in, t) 



e ^ 7, 7 — > e and 7 — > 7). Correspondingly the integral spectra have 



Bo 



dE n^{Eo,E,t) = 



Er. 



(17) 



These scaling properties of the approximation A solutions are a simple consequence of the absence of quantities with 
the dimension of energy in the shower equations. 
Energy conservation is reflected in the condition: 



dE E neiEo,E,t) 



dE E n^{Eo,E,t) = Eq 



(18) 



The solution of the shower equations for a monochromatic photon or electron [l[ is very simple for the the Mellin 
transforms of Ue and n^. The physically observable spectra can then be obtained inverting these transforms. The 
inversion can be easily performed numerically with a single path integration along a line in the complex plane, with 
exact results. Rossi and Greisen have also shown that it is possible to invert the trasform using a "saddle point 
approximation" , obtaining simple expressions for the spectra that are at the same time remarkably accurate and very 
instructive. 

The saddle point approximation solution for the differential spectra (valid for large t and E/Eq ^ 1) can be written 
as: 



naiEo,E,t) 



1 



1 



12-K 



with 



E 

lb'' 



E 
3t 



-(s+i) 



s=-s{E/Eo,t) 



t-2lii{E/E„) 

For the integral spectra the saddle point approximation solution is: 



Na{Eo, Ej-nin, t) 



1 



1 G„(s) 



V>W)t V Eo 



E„ 



(19) 



(20) 



(21) 



Note that both the electron and photon integral spectra diverge for i?min 0. For completeness a derivation of these 
well known results is sketched in appendix B, that also lists explicit expressions for the functions Ga(s). An example 
of the differential spectra in approximation A is shown in fig. [3l 

The saddle point solutions of the shower equations in approximation A exhibit several interesting properties: 



• For a fixed value (less than unity) of the ratio E/Eq (or Enun/Eo) the differential and integral spectra start at 
zero for t ~ 0, then grow with increasing t, reaching a maximum at the value tmax and then begin to decrease 
vanishing for t 00. The exponential factor e'^'i controls the t evolution of the solution. Therefore in good 
approximation shower maximum corresponds to Ai(s) 
finds the well known result: 



in 



E 
E'o 



(and therefore to s = 1). Using equation (^0)) one 

(22) 



In 



Eq 

E 



On can therefore rewrite equation (|20p in the form: 

E 



E< 



3t 



(23) 
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More in general, from the fact that the factor e'^^^'*^* controls the t dependence of the solution, it follows that 
to a good approximation one has: 



1 dNa 



Ai 



En 



or 



1 ^ ^ 

Ua dt 



s — ,i 

So 



(24) 



(25) 



The quantity s{E/Eo,t) is related to the shape of the energy spectrum around E, and has manifestly the 
meaning of the "local slope" of the spectrum at energy E: 



E 



dUa 

~dE 



E 
Eq 



(26) 



In other words the spectrum around E is well approximated as a power law E^^^^^\ Globally the spectrum 
is not a simple power law and the "local slope" changes as a function of E/Eq and t. For a fixed depth t the 
local slope s grows monotonically with E/Eq, as the spectrum becomes progressively steeper. For a fixed value 
of E/Eq the spectrum also becomes steeper with increasing t and as the local slope grows monotonically. The 
shape of the spectrum around energy E takes the form cx E~^ when the spectrum for this energy reaches the 



maximum. 



• Photons and electrons have spectra with very similar but not identical shapes, and accordingly the 7/e ratio 
changes slowly with energy. From the expressions for the functions Ga{s) (given in (jB14l - [Br7|) ) one finds: 



n^iEo,E, t) _ C(s) 
ne{Eo,E,t) (To + Ai(s) 



is) 



(27) 



That is the "asymptotic ratio" for the elementary power law solution with slope s — s{E / EQ,t). Note the 
remarkable fact that this result is independent from the nature of the particle (7 or e) that initiates the shower. 

• In fact the developments of showers initiated by a photon or an electron of the same energy are remarkably close 
to each other, demonstrating that the electron and photon population quickly tend to reach a sort of dynamic 
equilibrium feeding each other 

In approximation A there is no natural way to associate a single age to a shower because the equations do not 
contain any meaningful energy scale (except the energy of the primary particle). Therefore for a given depth t one 
can associate a different age to every ratio E/Eq (or i?min/^^o) according to equation ([20)1 . The age s describes at the 
same time the "stage" of the longitudinal evolution of spectrum (via equations ([24|) or ([25]) ). and the shape of the 
spectrum near energy E, that is well approximated by the power law cx E~^^^^\ The age also controls the 7/e ratio 
around E according to equation ((27| . 



III. ELECTROMAGNETIC SHOWERS IN APPROXIMATION B 



In approximation B the electron energy losses due to collisions are simply modeled as an energy independent loss e 
per unit of radiation length. The quantity e is the critical energy (in air e ~ 81 MeV). Accordingly, a term is added 
to the right hand side of equation ([5]) that describes the electron evolution: 



dnejE, t) 
dt 



dv ipo{v) 



ne{E,t) 



1 



E 



(28) 



^ du 



E 



t] +e 



dn^jE, t) 
dE 



The new system of equations [ (|28p and ([6])] does not have any more simple power law solutions of form ([7]). However, 
even in this case, one can introduce "elementary" solutions that have a constant shape in energy and evolve with t 
with the simple behavior e'^'^^*. Following Rossi and Greisen [l[ the elementary solutions can be written in the form: 



n^{E,t) = K e^^")* p(s,f) 
n^{E,t) = if e^(*)* ^-(^+1) 5(s,f) r^(s) 



(29) 
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that contain two additional functions p{s,x) and g{s,x). For large energy {E ^ e) the electron collision losses can 
be safely neglected, and the solutions coincide with the simple power law form of approximation A. This constraint 
tell us that the functions A(s) and r^{s) that appear in (j29p coincide with the functions discussed before and given in 
equations ([5]) and and that for large E/e the fmictions p{s,x) and g{s,x) asymptotically become unity: 

lim p{s,x) — 1, lim g{s,x) = 1. (30) 

X — ^oo X — yco 

Inserting expression (|29p in the shower equations one obtains two pairs of integro-differential equations (see appendix 
C) for the functions p{s,x) and g{s,x) corresponding to the two solutions for A(s). These equations can be solved 
numerically to obtain the functions Pi^2is: ^) a-nd gi^2{s, x). 

The physical meaning of the functions pi{s,x) and gi{s,x) is trasparent. If one injects power laws spectra of 
electrons and photons of form £'~(''+i) after a few lengths |A2(s)|^^ the spectra take asymptotically constant shapes 
given by: 



ne{E,t) = K e^i^-^)* pi (s,f) 



(31) 



(identical to but selecting the first of the two possible solutions) and continue their evolve in t as a simple 
exponential. 

The qualitative features of this asymptotic solution are easy to understand. The electron spectrum is a nearly 
perfect power law for 3> e but has a cutoff iox E ^ e, when electrons are absorbed because of ionization losses. 
The photon spectrum changes from a power law cx iJ^^'^+i) for i? ^ e to the form c>c i?^^ for i? <C e, reflecting the 
1/E dependence of the bremsstrahlung cross section. 

These physically intuitive properties are confirmed by the explicit calculation first performed by Rossi and Greisen, 
who have demonstrated [l| that the behavior of the functions p(s, x) and g{s, x) for a; ^ is: 

p{s,x) cx (32) 
g{s,x) cx x^ (33) 

The low energy behavior of the function p{s, x) implies that for E ^ Q the electron spectrum goes to a finite value. 
The energy integration of the electron spectrum therefore converges both for E oo (if s > 0) and for — > 0, and it 
becomes possible to talk about the total electron size. As expected the differential photon spectrum diverges cx E~^ 
for E ^ 0, and therefore the integral spectrum diverges logarithmically at the lower limit. 

The total electron size for the phenomenologically most important solution (that corresponds to Ai(s)) can be 
written as: 



/•CO I-CO / 

J dE ne{E) (X J dE E-^''+^'^ plis,— 

' dx pi(s, x) = e-' ^^li^iZ^ (34) 



The last equation defines the function Ki{s, — s). The physical significance of Ki{s, —s) can be understood comparing 
equation ([M]) with the integral 



dE E~^'+^'> 



£ 



s 



The electron size obtained integrating over all E the elementary solution (|29p differs from the integration of the simple 
form in the interval {e < E < oo) by a factor Ki{s,—s). Rossi and Greisen have shown how to calculate 

exact values of the function Ki{s, —s) for all integer values s > 0. For s = 1 one has Ki{l, —1) = 2.8948. A plot of 
the function A'l (s,—s) is shown in fig.IH 

The functions pi^2{s,x) and gi^2{s,x) can be calculated with numerical methods, as discussed in Appendix C. For 
a; > 1, it is also possible (see again appendix IC ip to express the functions as power expansion in 1/x with easily 
calculable coefficients. 

Figures 2] and [5] show (on a linear and log scale) the behavior of the electron spectrum plotted in the form 
due/dlnE — Ene{E) for three values of the index s. The integrated electron size Nf. is accounted for by parti- 
cles in the energy range 0.01 ^ E/e ^ 10. Note that an important limitation of the treatment in approximation B is 
that the electron mass is neglected. Accordingly the electron spectra extends down to E' to unphysical energy 
values below the electron mass. 
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A. Solutions for monochromatic electron or photon 



The solution of the shower equations in approximation B with the initial condition of a monochromatic electron or 
photon of energy Eq cannot be given with an exact closed form expression. Rossi and Greisen suggest to approximate 
the solution with the expressions: 



ne{-i)^e{Eo.,E,t) ~ [ne{j)^eiEo,E,t)]^X Pi 

n^(^^)_^{Eo,E,t) ~ [n^(^^)^jEQ,E,t)]^x gi 



E 

e 



(35) 
(36) 



with s(a;,t) given by (|20|) . These expressions combine the solution of the the shower equations in approximation A 
with the functions pi(s, x) and gi{s, x) introduced in section 2.2 as part of the "elementary solutions" to the shower 
equations. For E ^ e the solution coincides with the one obtained in approximation A, while for E ^ e the spectra 
have approximately the same shape of the "elementary solution" that corresponds to s{e/ Eo,t). 
In approximation B, it in possible and natural to consider the value 



E, 



3t 



t + 2 ln(£;o/£) 



(37) 



as the age of the shower. 

The total electron size of the shower Ne{EQ,t) obtained integrating over all energies is well approximated by the 
expression: 



7(e)- 



.(E'o, t) 



1 



Ea 
e 



Ki(s,-S) G'-y(e)^e(s) 



(38) 



- s=s(e/ Eo,t) 



where we have used the fact that integration over energy is dominated by i? ~ £ and the result . It can be easily 
seen, that the maximum of the size coincides with the condition Ai(s) — 0, that implies s — 1 and, solving equation 
(133: 



(39) 



E^ 

e 

Energy conservation in approximation B can be expressed with the equation: 



dE E ne{Eo,E, t) 



dE E n^{Eo,E,t) ^ Eo- e / dt' Ne{Eo,t') 



(40) 



The left-hand side of this equation is the energy contained in the shower particles at depth t, while the second term 
in the right-hand side gives the energy dispersed in the medium by the electrons as ionization. Equation (j40p also 
implies: 



/•oo 

/ dt N,{Eo,t) ^ Eo 
Jo 



(41) 



An example of the e and 7 spectra calculated in approximations A and B for a photon of initial energy 10"'^* eV at 
shower maximum is shown in fig. [6l The two solutions coincide for 3> e — 81 MeV, but deviate from each other at 
lower energy. The spectra in approximation B are strongly suppressed below the critical energy. The sum of the areas 
below the curves are proportional to the energy carried by each particle type. The curves of the approximation B 
solution enclose a smaller area because a part of the energy has been dispersed as ionization in the air (see equations 
(Ull) and (gOD). 

In approximation B showers generated by different primaries but having the same age s according to the definition 
([57]) have "essentially" equal spectra. This concept is illustrated in figure [71 that shows the electron spectra at 
shower maximum for showers generated by photons of different energy. The spectra are shown in two different 
representations. The first is of form {E dn/dE versus E), in this case the area below the curve is proportional to 
the electron multiplicity. In the other representation the spectra are shown in the form (i?^ dn/dE versus E), in this 
case the area below the curve is proportional to the amount of energy contained in electrons. In showers of the same 
age most of the particles have coincident spectral shapes, however the distributions of the highest energy particles 



9 



differ. High E particles account for a significant fraction of the energy contained in the shower, and are the reason 
why the evolution with t of a shower is not uniquely defined by the age, but depends also on the shower energy (or 
equivalently of the position t{s) where the age s is achieved. 

To summarize the results of this section: the explicit calculation of the average development of purely electromag- 
netic showers indicates that it is possible to define a shower age s: 



1 dN, 



dt 



3t 



3t 



t + 2ln{Eo/e) i + 2t„ 



(42) 



In showers of the same age the electrons and photons around and below the critical energy e (that dominate the total 
number of particles in the shower) have the spectra of the same shape. The energy spectra differ at larger energy. 



IV. THE GREISEN PROFILE 



The average longitudinal development of a purely electromagnetic shower generated by a photon or electron of 
energy Eq can be accurately described by a simple analytic expression introduced by Greisen f3| : 



^ Grciscn {Eo,t) 



0.31 



exp 



'('-2^°H^ + 2 1n(£;o/s) 



(43) 



The "Greisen Profile" is essentially identical to the more complex expression given in pT|) . The derivation Q of 
equation is simple and instructive, and requires the intelligent "recombination" of some the the results obtained 
above. 

The starting point of the derivation is the remark that the saddle point solution for the total electron size 
indicates the approximate validity of the relation: 



dN,{t) 
dt 



= Xiis)N,it) 



(44) 



with the age s given by 
equation (144]) as: 



One can now substitute for Ai(s) the approximation Ai(s) given in (|10p and rewrite 



dN.jt) 
dt 



Xiis)N,{t) 



3t 



t + 2tn 



1 - 3 log 



3t 



t + 2t 



The solution of this differential equation for the boundary condition N{t 

Ne{t) = ^max e~ 



max 



exp 



, 3, / 3t 
t ( 1 - - log 

2 y i 2 ^max 



N{t) . (45) 
is readily found as: 

(46) 



The normalization is fixed observing that the size at maximum for an electromagnetic shower can be obtained inserting 
the value s = 1 in equation ()38|) with the result: 



0.31 E^ 



where we have used: 



1 G 



"^^^^^ i^i(l,-l)= ^ 



27r yAf(T) 



2^ v/A^ 



^i(l,-l) 



0.31 



and 



(47) 



(48) 



"look" different from each 



Substituting this results in one obtains the final result (|43p . Expressions 
other but are essentially coincident numerically. 

A test of the accuracy of the Greisen profile solution can be performed verifying energy conservation using equation 
(|4T|) . This energy conservation condition is satisfied to better than 4.5% in the broad energy range from a few GeV 
to 1020 eV. 

In summary, the Greisen profile and the expression for the age in ([Ij (or 22) are equivalent to each other. The 
age definition IH) implies that the shower develops with the Greisen profile ([IS]) , viceversa the Greisen profile implies 
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the simple functional dependence for the age of equation ([T]) . The Greisen profile and the "Greisen age" ([T]) are (via 
the mapping Ai(s)) the integral and the derivative of each other: 

S(i,tmax) = -— r- A^(<) =Grciscn (t, imax) (49) 

A recent paper by Schiel and Ralston ^Tj] has complained that while the Greisen profile is shown and discussed 
many articles and textbooks, a full derivation is missing. The authors of make an attempt to "reverse-engineer" 
the steps performed orginally performed by Greisen to obtain this result, and arrive at the surprising and erroneous 
conclusion that the Greisen profile was "likely motivated by early numerical work in a time predating high-speed 
computer" . These comments miss the essential point that the Greisen profile is the result of the exact integration of 
a well defined (albeit product of some approximations) differential equation. 

The work [2] contains also a serious error when it argues that performing the derivative of the Greisen profile with 
respect to the critical energy and setting the value e ^ E, one obtains the electron differential spectrum at the energy 
E for a shower of primary energy Eq at the depth t. This error originates in a confusion between approximation A 
and approximation B for the shower equations, or perhaps more accurately in using the crude assumption to consider 
approximation B as nothing else than the introduction of a "sharp" cutoff for the electron spectra of approximation A 
for E <e. 

In fact performing the "trick" of the derivative with respect to s in the energy interval E ^ e. yield an interesting 
result that is proportional (but not equal) to the electron differential spectrum: 



rciscn 

{E^,e,t) 



de 



K{s,-s) neiEo,E,t) (50) 



(with s = s{t,Eo) according to equation ([H])). This result can be obtained comparing equations (PT|) and (|38p . When 
E aproaches e the interpretation of the derivative in l50l) ceases to be valid, and the electron spectrum takes the 
"universal" (age dependent) shape (x pi{s, E/e) £'~(*+i) with appropriate normalization. 

This last exercise is however instructive in the sense that it illustrates an important point. The age s determines the 
shape of the energy spectra of the "bulk" of the electrons and photons in the shower, but the distribution of the high 
energy particles depends on additional parameters. For the average development of an electromagnetic shower the ony 
additional parameter is the primary particle energy Eq (or equivalently the position of maximum tmax — lii(£'o/£))- 
The high energy particle content is crucial for the overall development of the shower. The Greisen profile implies 
that at each level t the shower contains a spectrum of high energy particles that is consistent with the shape of the 
development. This high energy particle content is not determined by the shower age, but (for each age) depends also 
on the primary particle energy Eq. The high energy {E ^ e) electron spectrum can be "extracted" from the shower 
profile via equation (|50p . 



V. UNIVERSALITY 

As discussed in the introduction, the recent works of Giller et al. [4] and Nerling et al. |5| have shown that for the 
same shower age (using the definition of equation ([1])), individual showers of hadronic primaries have electron spectra 
of the same shape. This result is clearly an important generalization of the result obtained in the previous section 
that the average development of purely electromagnetic showers. 

In the following we want to: 

1. Show that the shapes of the electron spectra calculated by montecarlo for individual hadronic showers in [1, Q 
are essentially identical to the Rossi-Greisen shapes, calculated for the average development of electromagnetic 
showers, and given for each age s by the expression pi{s, E/e) £'^(*+^'. 

2. Argue that this "universality" result is correct and expected, but that the definition of shower age should be 
modified from equation ([T]) that strictly speaking is only applicable to the average development of electromagnetic 
shower, to the much more general form ([3]). 

3. Show that also the shape of the photon energy distribution and its normalization relative to the spectrum of 
electrons in the shower are "universal" and determined by the shower age. 

As an example of the numerical coincidence of the Rossi-Greisen spectra of equation (|29p with the shape of the 
electron spectra calculated with montecarlo methods for the same value of the parameter s, in fig. [S]we compare the 
function pi(s, x) for the value s = 1 with the equivalent quantity (that is rie(s, E) £'(''+^)) from the fit to the electron 
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spectra at shower maximum obtained by Nerling et al. [5|. The two functions are nearly coincident, and agreement 
of comparable quality is obtained for all age values in the phenomenologically most important range 0.7 ^ s ^ 1.4 
(for more discussion see appendix |DJ . It is also interesting to note that the parametrization of Q for the electron 
spectrum has the form: 

\E + ai{s)\ [E + a2{s)\'' 

that has manifestly the same asymptotic behavior as the Rossi-Greisen shape at both low and high energy: 
ne{E) — > constant for ^ 0, and ne{E) cx £'"('*+i) for > e. 

The fact that the energy spectra calculated of individual hadronic showers coincide with remarkable accuracy with 
spectral shapes calculated for the average development of purely electromagnetic showers may appear at first sight 
surprising, but it is of course not a simple numerical coincidence, and has in fact a natural explanation. 

The first simple point is that in all shower types (electromagnetic, hadronic and also neutino induced) the total 
number of charged particles is essentially always (with the exception of very early and very late stages of development) 
dominated by electrons with energy around the critical energy e. In hadronic shower this happens because in each 
hadronic interaction a large fraction of the energy is transfered to photons via the production and decay of 7r° and r] 
mesons, these photons then generate the electromagnetic part of the shower that accounts for a growing fraction of 
the shower energy, and for most of the particles in the shower. 

The "universality" of the spectra in different showers of the same age, can then be immediately relating the "age" 
with the size slope A according to equation ([3]), and observing that the arguments outlined in the previous sections 
can be generalized to conclude that that the size slope A must be associated to a well defined shape for the "bulk" of 
the electrons and photons and to a well defined relative normalization between the two populations. 

For example, at shower maximum, when the shower size is "stationary" {dN/dt ~ 0), the photon and electron 
spectra must have spectral shapes and relative normalization that insure this stationarity of the shower size. A stable 
solution for this problem has been found in the previous section and is: 

[n,{E) cc E-^p{l,§) 

\ n^{E) cx 1.31 £;-2 5(l,f) ^ > 

To demonstrate formally that this in fact the general structure of the electron and photon spectra at shower maximum 
is in fact not trivial, however this is a very natural conclusion, observing how the result ([52|) emerge as the spectra 
at shower maximum for the average devolopment of the shower generated by both electrons or photons of arbitrary 
energy. 

Similarly, for each size slope A the quantity A^ (s) can be be identified with the parameters that labels the shape 
of the electron and photon spectra according to equation (|3ip . The argument that we have outlined to relate the age 
and the electron and photon energy spectra is independent from the nature of the primary particle, from its energy, 
and from the value of the depth where the shower is measured. 

The argument however does not allow to estimate the age from a closed form relation of type s(t,imax) such as 
equation ([1]). The fact that this is impossible can be easily illustrated with the example of the shower generated by a 
neutrino. In this case the ^max of the shower can of course be arbitrary large, and therefore for example the definition 
(fT|) returns for the entire shower development s ~ 1, that is of course meaningless, while it is physically transparent 
that the age concept mantain its validity and applicability, and the general definition ([3]) has no difficulty in dealing 
with neutrino induced showers. 

One may think that the neutrino example described above is "artificial" and that the problem that emerged can 
be "solved" for example shifting the origin of the depth measurement to the neutrino interaction point. However this 
point is in most cases unobservable, and this shifting procedure is operationally not well defined (and would open the 
problem of performing a similar shift of the origin of the t measurements also for hadrons and photon primaries). 

Another way of to see the problem for a closed form expression of the age of form s = s(i,t„iax) is that such a 
definition implies via equation (|44p the entire longitudinal profile of the shower. For example, as discussed before, 
the expression ([1]) implies that the shower develops with the Greisen profile (and viceversa). In general the 
longitudinal development of cosmic rays showers cannot be described accurately with the form and this failure 
implies limitations for the approximate definition p]). 



A. The "Gaisser Hillas" longitudinal profile 

The observations of longitudinal profile of high energy showers obtained with the detection of fiuorescence light, 
using the technique pioneered by the Fly's Eye detector, and currently in use by the HiRes and Auger collaboration. 
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are commonly fitted using a 4~parameters expression known as the "Gaisser-Hillas" profile: 



t-tn 



to 



exp 



to 



A 



(53) 



The maximum of this function is at t = tmax, where the size is equal to A^max, while to and A modify the shape. The 
Gaisser-Hillas profile (|53p implies the shower age: 



A 



(-1) 



1 



dNcait) 



Nonit) dt 



A 



(-1) 



1 {t ^max) 

'a [t - to) 



(54) 



where A^ ^\x) is the inverse function, of Ai(s). For completeness we note that the inverse of the function Ai(s) = A 

can has the explicit form: a[ ^''(A) — —3 ProductLog{— 1/3 exp[(— 1 — 2A)/3]}. To study the age near shower 
maximum, one expand in a power series around the position of the maximum, using as expansion parameter the 
quantity 5: 



^max) 



[tn 



to) 



(55) 



The first terms in the series expansion are: 



^-1+A- 



4A- 3 
4A2 



;A2 - 12A + 5 
8A=^ 



(56) 



For comparison, the age near shower maximum for the Greisen profile (|45p can be written as a power expansion in 
the quantity S': 



6' 



^max) 
^niax 



(57) 



with the result: 



St 



t + 2tr. 



oo 



k=l 



(58) 



A comparison of the expansions (|56p and ([55)1 shows explicitely how precisely expression ([55]) maps the true age of a 
shower. 



VI. CONCLUSIONS AND OUTLOOK 



The concept of the shower age can be very useful for the analysis of high energy cosmic ray data. The essence of 
the idea is very simple and can be summarized in a nutshell saying that the t-slope A and the i?-slope s of a shower 
are connected to each other by a one to one mapping. The i-slope (or size slope) is the fractional rate of change 
of the shower size with increasing depth (A = dN/dt). The i?-slope (or energy slope) is the integral slope of 
the (power law) energy spectra of photons or electrons above the critical energy. The mapping between A and s is 
given by A = Ai(s) ~ (s — 1 — 3 lns)/2. The spectra of photons and electrons have a more complex shape around 
and below e that is also determined by s (or A), and have a relative normalization also determined by s (or A). It is 
remarkable that the electron and photon spectra that correspond to different s (or A) have been calculated accurately 
with analytic methods by the pioneers Rossi and Greisen many decades ago. 

These properties of "universality" extend to the angular and lateral distributions of electrons and photons. This 
crucially important subject is not discussed here (see appendix lEl for some remarks). 

The definition of age discussed here is independent from the shape of the longitudinal development of a shower, 
and is therefore more general and accurate that the commonly used definition s ~ 3t/{t + 2tniax) that is correct only 
when the shower development is described by the "Greisen profile" . The average shape (and the fiuctuations around 
this average) of the longitudinal development of high energy cosmic ray showers is determined by the nature of the 
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primary particles and by the properties of hadronic interactions. The observation of these shape is a very important 
subject for future experimental studies. 

A definition of age that depends only on the derivative of the shower size can be applied also to neutrino-induced 
showers, and more generally is expected to remain valid for the description of all shower where the size is dominated by 
electrons. This includes showers generated by exotic primaries, or the presence of unexpected physics (or unexpected 
fluctuations) in the development of the showers by primary particles of known nature. The search for events that 
have unusual longitudinal developments, such as multiple maxima is an interesting direction of research. It is likely 
(and at least the best possible a priori assumption) that the spectra of the electromagnetic component around and 
below the critical energy will, also in these cases, be controled by the shower age. 

These ideas can be useful in the analysis of high energy cosmic ray observations in several ways. As examples: 
(i) the knowledge of the variations of the electron energy spectrum during the evolution of a shower can be used to 
obtain a better reconstruction of the longitudinal profile of the shower in observations that use fluorescence and/or 
Cherenkov light detectors (see for more discussion); (ii) the reconstruction of the shower age from the lateral 
distribution of its electromagnetic component can in principle help in the reconstruction of the energy in surface array 
measurements; (iii) in case of hybrid measurements of the showers, the redundant measurement of the age (from the 
size longitudinal development and the lateral distribution of the electromagnetic component at the ground) can allow 
to disentangle a muon component, allowing composition measurements, or test of hadronic interaction models (see 
[To| for more discussion). 

It should finally be stressed that the "universality" in the electromagnetic component of high energy showers, is 
clearly an important analysis tool, but gives only a partial information about the shower. Other information is 
contained in the shower muon component, moreover the shower core, that is essentially undetected in large area 
shower arrays, in most cases also contains a significant amount of energy. The energy contained in the core must be 
"infered" from the information obtained at large distances from the shower axis, introducing unavoidably some model 
dependence in the energy reconstruction. 
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APPENDIX A: CROSS SECTIONS FOR FUNDAMENTAL PROCESSES 



In this appendix we list the expressions of the differential cross sections for the bremsstrahlung and pair production, 
the fundamental processes that control the development of electromagnetic showers. It is convenient to measure the 
column density X in units of radiation length Xq (with the notation t = X/Xq). The radiation length in air [T^ is 
approximately 36.66 (g cm^)^^. 

The probability per unit of radiation length that an electron of energy emits a photon of energy E^ = v E^, has 
the asymptotic form: 



1 



ip{v) = - 



1- ( --26) (1- «) + (!- «r 



(Al) 



For the pair production process 7 — > e e^, the energy distribution of the final state electron is: 

V'(u) = (1 - uf + ( ^ - 26 ) (l-u) u + u^ 



(A2) 



with u = E^-/E^. Integrating over all u values one obtains the pair production probability per radiation length: 



7 



9 3 

In the previous equations 6 depends on the atomic number of the medium: 

1 



18 log(183Z-i/3) 



(A3) 



(A4) 



For air one has 6 ~ 0.0135. 

Important momenta of the functions ip{v) and ip{u) are: 



A{s) 



dv ip{v) [1 - {1-vY] 



26 



r^(i + -'^) 
r(i + s) 



+ 7 + 



s ( 7 + 55 + 126 (2 + s)) 
6(1 + s) (2 + s) 



B{s) 



Cis) 



du u'^ 'ipiu) 



dv w'* fiv) 



2 (l4+lls + 3s2-66 (1 + s)) 
" 3 (1 + s) (2 + s) (3 + s) 

8 + 7s + 3s2 + 66(2 + s) 



3s (2 + 3s + s2) 

In equation (jA5[) r'(z)/r(z) is the digamma function and 7 is the Euler gamma constant 7 ~ 0.577216. 



(A5) 
(A6) 
(A7) 



APPENDIX B: SHOWER EQUATIONS IN APPROXIMATION A 

In this section we sketch a derivation of the solutions of the shower equations in approximation A for the initial 
condition of a monochromatic electron or photon. 

The first step is to introduce the Mellin transforms Me{s, t) and M^{s, t) of the electron and photon spectra ne{E, t) 
and n.y{E, t). In general the Mellin transform of the function f{E) is defined as: 

poo 

M/(s)= / dEE' f{E) (Bl) 

with s a complex parameter. The Mellin transform converges in a strip bounded by two straight lines parallel to the 
imaginary axis (si < 3fJ[s] < S2). The inverse transformation is: 

f{E) = I ds M}{s) (B2) 
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where the integration path C runs parallel to the imaginary axis within the strip of convergence of Mf(s). 
Applying the operator: 

dE E" 





to the shower equations ([5]) and® one obtains a system of two linear differential equations for Me(s, t) and M^{s, t): 

^M^hll ^ ~A{s)M,{s,t)+B{s)M^{s,t) (B3) 

+C(s) Me(s, t) - (To M^(s, t) (B4) 



dM^{s,t) 



dt 

The general solution of this system can be easily obtained as a linear combinations of the exponential functions e^^'^^^ * 
and e^2(")*. 

The shower generated by an electron of energy Eq at depth < = corresponds to the initial condition: 



n,{E,Q) = 5[E~Eo\ 
n^{E,{)) = 



Me(s,0) = [E^y 
M^(s,0) = 



while the shower generated by an initial photon of energy Eq corresponds to the initial conditions: 

ne{E,Q) ^ 
n^{E,{)) = 6[E-Eo\ 



or: 



Me(s,0) = 
M^(s,0) = {EoY 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



Using these boundary conditions one finds the solutions: 

M (En 'it) - —£^^LEo_fp>^l{s)t_ \2{s)t^ 

iVJe— — Ai(s)-A2(s) e J- 

and 

m^_,e[-C^O,b,l) - p^^^ Ai(s)-A2(s) e I 

The functions ne,-y{EQ, E,t) can be obtained inverting the Mellin transformation using (jB2[) . or more explicitely: 

1 /• 1 pSQ+i CO 

n^{E) = — ds M„(s) = — ds M„(s) (Bll) 

27^1 Jc 27ri 

(with the subscript a than runs over the 4 cases: (e — > e), (e — > 7), (7 — > e) and (7 — > 7)). This integral cannot be 
done exactly analytically, however with modern tools it is trivial to obtain the numerical result with any desired level 
of accuracy. In fact the integrand function is well defined (and available in computer libraries) for all complex values 
s. The imaginary part of the integral vanishes while the real part gives the physically observable spectra. 

Rossi and Greisen have also shown that using the "saddle point" approximation it is possible to obtain simple 
analytic expressions that are a very good approximation of the exact results for t not too small, and that remain very 
instructive and useful. The basic idea behind the saddle point approximation is that the integrand in (jBlip is an 
analytic function in the variable s. For any analytic function f(z) = f{x + iy) one has: 

— - H = 
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This implies that if the function /(z) has a minimum for a real value z when z runs along the real axis, the function will 
then have a maximum at the same point along a path that is at right angle with respect to the real axis. The integral 
is then dominated by the value of the function near the maximum, and can be performed analytically approximating 
the integrand with a Gaussian function, and using the well known result that if Q{x) is a quadratic form: 



1 2 
91 a; + -92 a; 



with coefficient 52 > 0, one has: 



+00 



exp {-Q{x)\ 



(B12) 



(where x = —qxjq^ is the point where the quadratic form Q{x) is minimum, and the integrand is therefore maximum). 

One can now apply this idea to the integral (|B11[) . For t not too small, one can neglect the term proportional to 
exp[A2(s) i] in the expressions for the McUin transforms given in equations (|B9p and (|B10[) and rewrite the integrand 
of the inverse transform as: 



27r« 



where the functions Gq(s) are: 



2TTi E 



E 
Ea 



G^^eis) 



G 



1^1 



(^) 



[gp + Ai(s)] 
Ai(s)- A2(s) 

G{s) 
Ai(s)- A2(s) 

l_ [gp + Ai(5)] [o-Q + A2(s)] 

G{s) Ai(s)-A2(s) 

kp + A2(5)] 
Ai(s)~ A2(s) 



(B13) 

(B14) 
(B15) 
(B16) 
(B17) 



In equation (|B13|1 the integrand of the inverse Mellin transform has been written as the product of a function that 
changes rapidly with s (in square parenthesis) and a function that is considered as slowly varying. The part of the 
function that is rapidly varying with s has a minimum along the real axis for the value s determined by the implicit 
equation: 



d_ 

ds 



E 



,Ai(s)t 



= 



or equivalently: 



A'(s)i + ln 



= 



(B18) 



(B19) 



This equation has an explicit solution if one substitutes for Ai(s) the Greisen analytic approximation Ai(s) given in 
(fTOl). The solution is: 



E 

:t 

SP 



it 



t - 2 \iy{E/Eq) 



(B20) 



One can now complete the calculation at the saddle point s — s approximating the integrand as a gaussian function 
and using equation (jB12[) with 92 — A"(s) t with the result: 



1 1 

na[Eo,E,t) ~ — ^= 
^p v27r 



E 



(B21) 



J s=s(E/Eo,t) 
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The integral distributions can be calculated noting that for t not too small the integration is dominated by E close 
to i^niin- Neglecting the slow variation of s with E one finds: 



-^a(£^inin, Eq, t) 



1 



1 G„(s) 
s 



(B22) 



J s=s(£;„,i„/£;o.t) 



An alternative method to obtain equation (|B22p is to observe that if is the integral of the function f{E) 

for E > i^niini then the Mellin transform Mp{s) is given by: 



MF{s) = -^Mf{s + l) 
and performing the inversion with the the saddle point approximation. 



(B23) 



For completeness we note that Rossi and Greisen suggest some slightly more complex forms for the saddle point 
solutions in place of equations (|B21[) and (|B22[) for a better approximation with the exact solution. The idea is to 
make a better choice for the "fast varying" part of the Mellin transform. This can be done introducing the quantity 
TOq and rewriting the decomposition (|B13[) as: 



~" Zm E 



2711 



E 



(B24) 



The esponent rria is chosen "ad-hoc" for better quantitative results. Rossi and Greisen [l| suggest the values: 

(B25) 
= A (B26) 



''e — *7 



'TO 



7 — >e 



The part of the function that is considered as varying rapidly with s is indicated in square parenthesis in (|B24p and 
has a minimum along the real axis for the value s determined by the implicit equation: 



E: 



A'(s)t + ln^ +^=0, 



(B27) 



that can be solved explicitely if one substitutes Ai(s) with the Greisen analytic approximation. The solution is: 

3t ^ 2 rria 



t~2 \n{E/E^) 



(B28) 



One can proceed with the saddle point solution, noting that for the Gaussian approximation of the rapidly varying 
function the parameter q2 is now given by: q2 ~ A"(s) t — TOq/s^. The differential spectra can then be written as: 



^0 v2tt 



y/X"{s)t-ma/s^ 



E 
E'o 



-is+l) 



,Xl{s)t 



(B29) 



J s=Sa(E/Eo,t) 



Similarly a more complex expression can be written to improve on equation (|B22p for the integral spectra. For large 
Eq and large t one can neglect the terms nia/s^ and the more complex expressions for the differential and integral 
spectra coincide with the simpler results (jB2ip and (|B22p . 



APPENDIX C: ELEMENTARY SOLUTIONS IN APPROXIMATION B 



The form of the solution for the shower equation in approximation B is given in equation (|29p . Inserting this 
expression in the shower equations one obtains two integro-differential equations for the functions p{s,x) and g{s,x): 



A(s) pis,x) 



2C{s) 
era + A(s) Jq 



du u^'il){u) g ^s, — ^ 
p{s,x) - (1 - vfp ( s 



+ 



dv ip{v) 

p{s,x) dp{s,x) 



(CI) 



dx 
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and 



is,.)^^j\vv^p{s,^). (C2) 



Cis) 

It is possible to eliminate the function g in (jCip . obtaining an equation for p(s, x): 



\{s)p{s,x) = TVTT / duu^'ipiu) [ dvv" (p{v) p (s, —) 



dv (p{v) 







p(s,a;) - (1 - vYp s, 



X 



1-v 



(C3) 



The two solutions for p{s,x) are obtained substituting in (|C3[) the two solutions for A(s) given in ([H]). Note that 
equation (|C2p is a simple integral, therefore from the knowledge of the electron energy spectrum, that is the functions 
Pi, 2(5, x), it is trivial to obtain numerically the corresponding photon spectrum, while it is not entirely trivial to solve 
numerically equation (jC3p for p{s, x). 

1. Power Expansion for the functions p{s,x) and g{s,x) 

A useful result obtained by Rossi and Greisen is the demonstration that the functions p{s,x) and g{s,x) can be 
expressed as a power series in 1/x: 

p{s,x) = ^c,(.)x-^- = l + ^ + H^ + ... (C4) 
g{s,x) ^ yd,{s)x-^ = l + ^ + ^ + ... (C5) 



with easily calculable coefficients. In fact, inserting these power series forms in equations (|C3|) and (|C2p one obtains 
simple recursive relations for the coefficients: 

c„_i (s + n) 



where the function -F(s, n) is: 



^, ^ ^ / N ^/ X B{s + n)C{s + n) 

F(s,n) = A(s)+A(s + n)- ^ ^-^ (C8) 



and the functions A(s), B{s) and C(s) are listed in appendix Kl Clearly a power series in 1/x solution cannot be 
used for x close or below unity (that is for E close or below the critical energy e), however this expansion allows to 
test in its range of validity solutions obtained with different methods. 



APPENDIX D: ANALYTIC AND MONTECARLO SOLUTIONS 

The calculation of the spectra of electrons and photons in shower of the same age obtained by Rossi and Greisen in [l[ 
and discussed in this paper has the limitations that are intrinsic to the realistic but simplified theoretical framework 
(approximation B) that has been used. This framework introduces several simplifications: the cross sections for 
bremsstrahlung and pair production have always the asymptotic form that is strictly speaking only valid at very high 
energy, the electron collision losses are treated as a simple energy independent constant, and Compton scattering is 
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entirely neglected. Also in approximation B the electron mass is neglected and the electron spectrum extends down 
to zero energy. A montecarlo calculation of the spectral shapes can of course avoid all these limitations and is in 
principle more accurate, even if it has its own limitations and difficulties 

A comparison of the Rossi-Greisen shapes with the numerical results of [1, [E| > shows remarkable agreement but 
also some small differences, that could be interesting to explore in more detail. As an illustration, normalizing the 
high energy spectra to ne{E, s) E^^'^^^\ the quantity 

/•oo 

s / dEne{s,E) (Dl) 

is given by the function Ki{s, —s) for the Rossi-Greisen calculation. The numerical integration of the Nerling et al 
parametrization 5] gives results that differ by 5-10% (at shower maximum the difference is 4.5%). A comparison of 
the results is shown in fig. [9l 

As an additional test we have calculated the first coefficient in the development: 

pi =E'+^n,{s,E)^l + c^{s)^+C2{s) (|)' + ... (D2) 

For the Nerling parametrization ([51]) the first coefficient is ci(s) = — [ai(s) + sa2(s)], while for the Rossi-Greisen 
solution the coefficient is given by (|C6[) . A comparison of the two estimate is shown in fig. (TU] 

The origin of the (small) differences between the analytic and montecarlo solutions merits further studies. A 
possible explanation is a more precise description of the physics of the electromagnetic interactions in the montecarlo 
calculation. 



APPENDIX E: LATERAL DISTRIBUTION 



It is intuitive that in showers of the same age "most" electrons and photons have not only the same energy 
distributions but also the same angular distributions and, for the same density profile of the medium where the 
showers are propagating, also essentially equal lateral distributions around the shower axis. 

The problem of calculating the electron lateral distribution has attracted considerable in the past. Nishimura 
and Kamata [3] solved numerically the 3-dimensional shower equations in approximation B to obtain the (energy 
integrated) lateral distribution of electrons propagating in a medium of constant density. Their result were fitted by 
Greisen [s^] with the approximate form: 

p,[s, r) = TVe K{s) \ f-) f 1 + (El) 

with rg is the Moliere radius and K{s) = (27r)~^ r[4.5 — s]/(r[4.5 — 2 s] r[s]) is a normalization factor. After these 
works several other authors have given different parametrizations of the lateral distribution as a function of an age 
parameter (see for example Hillas in [ll| . 

In the view of this author it is in fact not possible to have a single parametrization for lateral distribution, because 
it is essential to consider the properties of the detector that measure the shower at the ground. For example the 
results of Nishimura and Kamata refer to the total number of electrons integrated down to zero energy, however in 
most cases the detectors of shower arrays do not sample the electron number but an energy deposition, and in any 
case it is always needed to take into account some contribution from photons in the shower. Therefore one needs to 
combine appropriately the electron and photon contributions with the detector response. An additional complications 
is of course that for hadronic primaries one has to disentangle the the electromagnetic and muon components of the 
shower. 
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FIG. 2: Equilibrium photon/electron ratio: ri^\s) = C{s)/[ao + Xi{s 
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FIG. 3: Electron and photon energy spectra calculated in approximation A for the shower generated by a primary photon of energy 
10^® eV at three values of the depth (t = 14.6, 23.7, t=41.3 that correspond approximately to age s = 0.7, 1 and 1.4). Thick (thin) lines 
are for electrons (photons). The spectra are shown in the form "(^) versus E. The area below each curve is proportional to the 
amount of energy transported by each particle type at the depth considered. 
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FIG. 4: Plot of the energy distributions of electrons calculated in approximation B for three values of the age parameter s (s = 0.7, s = 1 
and s = 1.3). The distributions (in the form Erie = drie/dlnE) are calculated as: pi{s,E) and are renormalized to have a total size 
of one electron. 




FIG. 5: Plot of the energy distributions of electrons and photons calculated in approximation B for three values of the age parameter s 
{s = 0.7, 1 and 1.3) The distributions (in the form Eriej = drie^-y /din E) are calculated as: p(s,E) E^"^ for electrons g{s, E) E~''rj(s) 
for photons and renormalized to have a total size of one electron (and the correct 7/e ratio) 
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FIG. 6: Comparison of the approximation A and approximation B solutions for the electron and photon spectra in the shower generated 
by a primary photon of energy 10^* eV at shower maximum (t = 23.7). The area below each curve is proportional to the amount of energy 
transported by each particle. 
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FIG. 7: Electron spectra for the showers generated by primary photons of energy 10-'^'' cV, 10^** eV, and 10^" eV at shower maxinmm 
{t = 18.6, t = 2.3.7 and t = 27.8). The normalization is chosen so that the spectra are equal at i? = £ = 81 MeV. The top panel shows 
the spectra in the form E n{E) versus E; the bottom panel shows the same spectra in the form E^ i^{E) versus E. The similarities and 
differences between the curves illustrates the concept and the limitations of the "universality" of the electron spectra in showers of the 
same age (in this case s = 1). 
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FIG. 8: Plots of the functions p{s,x) and at shower maximum (s = 1). The points show the results obtained taking the first 

4 terms the power series developments l|C4[l and (|C5| . The dashed line shows the corresponding fit of Nerling et al 5] for the electron 
spectrum at shower maximum. 




2.0 



FIG. 9: Plot of the function Ki{s, —s). The dashed line shows the normalization of the electron spectrum of Nerling et al 




FIG. 10: Plot of the coefficient ci(s) of the expansion ||C4[ | for the pi{s,x) function. The red dashed curve is the coefficient for the 
numerical calculation of Nerling et al. 



